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Abstract 

The time-dependent behavior of a two-level system interacting with a quantum oscillator system 
is analyzed in the case of a coupling larger than both the energy separation between the two levels 
and the energy of quantum oscillator (f2 < oj < A, where f2 is the frequency of the transition 
between the two levels, u is the frequency of the oscillator, and A is the coupling between the two- 
level system and the oscillator). Our calculations show that the amplitude of the expectation value 
of the oscillator coordinate decreases as the two-level system undergoes the transition from one 
level to the other, while the transfer probability between the levels is staircase-like. This behavior 
is explained by the interplay between the adiabatic and the non-adiabatic regimes encountered 
during the dynamics with the system acting as a quantum counterpart of the Landau-Zener model. 
The transition between the two levels occurs as long as the expectation value of the oscillator 
coordinate is driven close to zero. On the contrary, if the initial conditions are set such that the 
expectation values of the oscillator coordinate are far from zero, the system will remain locked on 
one level. 

PACS numbers: 42.50.Md, 42.50.Hz,63.20.Kr, 85.25.Cp, 85.85.+j 
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One of the most studied quantum-mechanical models is the two-level system interacting 
with a quantum oscillator. It is used in a wide range of phenomena, especially in atomic 
physics where it describes a two-level atom coupled to a quantized electromagnetic fields 
A challenging enterprise is to extend the model to "artificial atoms" in a condensed-matter 
environment. These small solid-state devices like flux lines threading a superconducting 
loop, charges in Cooper pair boxes, and single-electron spins exhibit quantum-mechanical 
properties which can be manipulated by currents and voltages.—^ 

The solid-state devices offer wider regimes for the coupling strength between the two-level 
system and the oscillator. Typically, the coupling strength in atomic systems is X/uj = 1CT 7 — 
10~ 6 . 2 Similar dipolar coupling in Cooper pair boxes and Josephson charge qubits is 3 — 4 
order of magnitude larger than the coupling in atomic systems.—^ In contrast to the dipolar 
coupling, the capacitive and inductive couplings show even larger coupling strengths . 8 ' 9 ' 10 
Recently it has been argued that values X/uo « 1 are possible to be achieved experiment ally- 
In this brief report we examine the regime X/uj > 1. We will show that, although the 



oscil 



ator dynamics "follows" the dynamics of the two- level system as in the case studied in 



Ref. ljj the general features are different: the system undergoes the transition from one level 
to the other with a sudden change in the transition probability, whenever the expectation 
value of the oscillator coordinate is close to zero. 

The Hamiltonian of the system is written as (K — 1): 

v 2 + oj 2 a 2 

H = 2 +\qcr z + na x , (1) 

where a z , and a x are the spin operator matrices, and p and q are the oscillator coordinates. 
The essential parameters are u, X, and Q, associated with the frequency of the oscillator, the 
coupling strength of the two-level system with the oscillator, and the splitting frequency of 
two-level system, respectively (our parameter A is scaled up by a factor of 2\/2 with respect 

11). In Ref . Q it was assumed that the splitting frequency is much 



to parameter A in Ref. 



smaller than the frequency of the oscillator, thus the problem can be cast into the displaced 
oscillator basis which is the basis for the first two terms of the Hamiltonian (0). This 
displaced oscillator basis is found by applying the unitary transformation U = exp (^jrcr 2 ) 
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to the Hamiltonian (JIJ). The transformed Hamiltonian become; 



H' = + % [a+exp (-i*J) + - |S > ( 2 ) 

with o- + ((T~) as the spin-^ creation (annihilation) operator, than time- dependent or/and 
time-independent perturbation calculations can be performed as long as Q/u « 1. A 
similar path was followed by Schweber.— Perturbation calculations on the energy spectrum 
have been performed and compared with the exact calculations in Fig. 3 of Ref. The 
authors extended the comparison to the regime VL/uj > 1 and they found some quantitative 
and qualitative resemblance between the exact solution and their approximate solution (Figs. 
3(b) and 3(d) in their paper). The quantitative and qualitative agreement shown in Fig. 3 
of Ref. Q can be explained in simple terms as follows. The net effect of exp (Ipf) on the 
wave function is to displace it by the amount jpy. Thus, the effective splitting given by the 
second term of the Hamiltonian (J2J) will be quenched, such that perturbation calculations 
on the Hamiltonian (|2j) can be extended to larger values of Q. 

Irish and coworkers 11 studied the collapse and revival of the wave function for Vt/uo « 1 
and X/u < 1. In the following we will explore the dynamics of the very strong coupling 
regime {X/uo > 1 and Vt/u < 1). The Hamiltonian (J2J) would be able, in principle, to 
explain the dynamics in the very strong regime because the effective coupling between the 
two wells generated by the displaced oscillators is quenched by the separation of these two 
potential wells. The effective coupling between the two wells decreases exponentially^ with 
the strength of A as it can be shown also by analyzing Eq. (J2J). However, to gain a new 
insight into the dynamics of the very strong regime, we follow a different approach. We 
perform another unitary transformation , 12 ^ U = exp (iA (q) a y ), on the Hamiltonian (flj). 
with tan (A (q)) = — j^, to obtain adiabatic motions which are valid for either very strong 
coupling A >> Q,u (case A) or large level splitting Q >> A, a; (case B). The intuitive 
picture of the transformation is a g-dependent rotation around the y-axis that brings the 
effective field seen by the two-level system along the z-axis. The Hamiltonian resulting from 
the above unitary transformation is 

h> = W + uV) + §^Sw + *- 

The above Hamiltonian Q will be used below for our investigation. The last term (cr^-term) 
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FIG. 1: The shape of the two sheets of the adiabatic potential generated by the Hamiltonian in 
Eq. (flf). The solid line is the lower sheet and the dotted line is the upper sheet. The second term 
in Eq. (jlj) is sharply peaked around origin and it is plotted with dashed line. 

in Eq. (j^J is small as long as Xq >> Q (that is supposed to be fulfilled in Case A) or Xq << Q 
(that is fulfilled in Case B). The inequalities Xq » Q and Xq << Q should be understood 
as operator inequalities in the sense that they have to be satisfied as inequalities for matrix 
elements in a certain basis. The cr y -term is the non-diagonal term and it accounts for the 
non-adiabaticity. Without the <7 y -term, the Hamiltonian Q reveals an adiabatic motion 
with one part (either the two- level system or the oscillator) becoming fast, while the other 
part becoming slow. 1 - Thus, the adiabatic motion is generated by the first three terms, 



In the very strong coupling case (A >> ft, u), there are two adiabatic potential sheets 
coupled by the a y term . 12 ' 14 The shape of the adiabatic potentials is presented in Fig. [TJ 
Dynamically, the oscillator is fast and the two-level system is slow.— The lower adiabatic 
sheet has two minima located at 



Q 2 X 2 



(4) 
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with the value 



V (q min ) 



M 2 

2X 2 



8w 2 ' 



(6) 



4 



exact 

adiabatic 

oscillator 






1.0- 




y(2m)= -5,0/(0=0.1 


A 








D N 
V 


0.5- 


(b) 




Ot 








0.0- 









15-j 




10: 


A 




cr 
V 


5: 




0: 




1.0- 


A 




b N 


0.5^ 


V 


CM 






0.0- 



adiabatic 

oscillator 




10 20 30 40 50 60 70 80 

cot 





V(2ra)= -5, D 


0=0.5 


1 1 1 1 1 1 


— ■ — 1 — ■ — 1 — • — 1 — ' 





FIG. 2: Time-dependent behavior of the expectation value of: (a) q— coordinate, with — = 0.1 
and = —5; (b) 2a z , with — = 0.1 and = —5; (c) q— coordinate, with 2 = 0.5 and = —5; 
(d) 2a z , with ~ = 0.5 and ^ = —5. The solid line denotes the exact dynamics, the dashed line 
represents the adiabatic dynamics, and the dotted line shows the corresponding displaced oscillator. 
For convenience, we have chosen negative values of A. We notice that the strength of the coupling 
is, actually, ^. 



and its second derivative at the minimum points 



u* 1 - 



(7) 



The lower adiabatic sheet is very close to the unperturbed (Q = 0) displaced harmonic 
potential sheets as long as A >> Q, u. Therefore one might expect similar dynamic behavior 
as the one studied in Ref. Q. However, the dynamics is different for very strong coupling. 
In Fig. |21 we show the dynamics of the expectation values of g-coordinate and a z . The 
details of numerical integration are given in Ref. We compare the exact dynamics with 
the dynamics of the adiabatic potential (Eq. and with the dynamics of a displaced 

oscillator (Q = 0). One can notice in Fig. [2] that the amplitude of the expectation value 
of g-coordinate decreases with time. It occurs as the system undergoes the transition from 
one level to another. In the same time, the transfer probability is staircase-like. Moreover, 
the slope of the amplitude decrease depends on Q but the frequency of the oscillations does 
not change significantly. Although it is weaker, the same quenching of g-coordinate appears 
for the adiabatic motion generated by Eq. (JU). It is weakly dependent on Q and it occurs 
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as the system undergoes the transition from one potential well to the other potential well. 
This was pointed out by Wagner— who showed that the transition rate from one well to the 
other is dependent on Q in second order in the adiabatic approximation. In contrast to his 
paper 12 , Fig. |2] shows clearly that the adiabatic motion is not an accurate description of the 
full dynamics. This behavior will be explained below as interplay between the adiabatic and 
the non-adiabatic regimes encountered during the dynamics. 

In order to explain the dynamics, we employ the equation of motion in the Heisenberg 
picture for an operator, = i [H, A]. Taking the expectation value, the equation becomes 
4:{A) = i([H,A\). We apply this last equation to the Hamiltonian represented by Eq. ((TJ). 
The corresponding equations are 

i<?> = <p>> 

f t (p) = -^( q }-X(a z } 1 

f t (cr x ) = -\(q<T v ), 
j-A(jy) = X(qa x ) -Q(<T Z ), 



Eq. (jHJ) is an infinite chain of coupled ordinary differential equations. The chain can be 
broken by making assumptions like 

fay) = (q){<7y), 
(qa x ) S (q)(a x ). 

The approximations made in Eq. Q are valid if one part (the oscillator) is fast and the 
other part (the two-level system) is slow as in the usual adiabatic approximation . 15 > 16 Thus, 
combining (jHJ) and Q, one can show that Eq. (jHJ) can be approximated and then recast as 

Hi) = <p>, 

i{p) = -u\q)-\{\h\ 2 -\h\ 2 l (1Q) 

with fi (/ 2 ) as being the probability function of the level 1(2). In other words, |/i| 2 (|/2| 2 ) 
is the probability of the system to be on the level 1(2). 
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FIG. 3: (a) Comparison between the dynamics of (cr z ) (solid line) and \{q)/Cl (dotted line). It 
shows that (a z ) changes whenever \(q)/£l is close to 0. The upper panel shows the scale around 
origin, while the lower panel shows the full scale of the plots, (b) Pictorial explanation of the 
dynamics of the two-level system in this very strong coupling case. The transition from the left 
well to the right well occurs at the origin where the adiabatic condition (\{q) >> O) is not fulfilled. 

Eq. (fTU|) sheds a better light on the relationship between the oscillator and the two- 
level system. The first two equations in (fTU|) are the equations of the classical harmonic 
oscillator displaced by the amount §(|/i| 2 — |/2| 2 )- This will imply that the time variation of 
(|/i| 2 — I/2I 2 ) = 2(cr z ) modulates the amplitude of the oscillator as it can be seen in Fig. [21 
The last two equations in (fTU|) explain the time-dependent behavior of (|/i| 2 — l/^l 2 )- First, 
let us assume that X(q) >> Q, which is the adiabatic condition and it is supposed to be true 
most of the time for A >> Q, u. The last two equations in (fTU|) will be approximated by 
i f l = I (g) f l and % 4 f% = — ~ (q) ^2- This means that the probability functions fi and 
f2 acquire just a phase factor. Therefore, there is no mixing between fi and f2, and (<r z ) is 
constant in time. Now, let us assume that the reverse is true, X(q) < Q. Then, the last two 
equations in (fTU|) will be approximated by i 4 fi = 2-/2 and i f2 — he., fi and will 
mix, and (a z ) will change. These assertions are proven numerically in Fig. Ela), where one 
can see that (a z ) changes whenever \{q)/Q = 0. 

In Fig. Efb), we give a simple explanation of the dynamics of the transition from level 1 
to level 2. Suppose that the wave packet of the system is in the left well which is associated 
with level 1. As soon as the wave packet reaches regions close to the origin it "sheds" a 
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FIG. 4: (Color online) Dynamics of (a z ) with different initial conditions. The arrows indicate 
various initial conditions in terms of the displaced ground state wave-function of the initial oscil- 
lator. The oscillator displacements are — 0.5(/o (it produces the faster transition), 0, 0.5(/o> l-5go, 
and 2.5go (it locks the system on level 1), with go = £j = 0.2 and = 3. 

part of itself into the right well and the rest of it returns into the left well. Basically, the 
transition from level 1 to level 2 occurs during short periods of time when the system is 
out of its adiabatic regime. The picture presented here is closely related to Landau-Zener 
theor y 1 ^ 1 ? that treats a quantum two-level system placed in a slowly varying external field. 
Near the crossing point the adiabaticity is violated and the system can escape from the state 
it occupied initially to another one. In the present case, the coupling is purely quantum 
through the quantum oscillator, hence the oscillator explores all possible trajectories (in the 
sense of Feynman's path integral) and the ones that explore the region at the level crossing 
can induce a non-zero nipping probability for the two-level system. 

In Fig.|3]we show the dynamics of (a z ) as depending on initial conditions. We consider as 
initial conditions various displaced ground state wave-functions of the harmonic oscillator. 
The displacements are given in terms of go = — y^, WIUC h 1S the displacement of the unper- 
turbed oscillator (Q = 0) corresponding to the level 1. Assuming that the system is initially 
on the level 1, its g-coordinate, (q), will tend to oscillate around go, starting from the initial 
displacement. Fig. 0] shows that the system undergoes faster transitions from level 1 to level 
2 for those initial conditions which generate dynamics that fulfills the condition X(q) << Q 
for longer time. In addition to that, the system can be locked on level 1 as long as the initial 
expectation value of q is set to be close to go. 
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In conclusion, we investigated the dynamics of the two-level system interacting with a 
quantum harmonic oscillator in the regime of very strong coupling (A > to > Q)). This 
regime generate an adiabatic motion defined by the condition X(q) » Q. We have found 
that the amplitude of the expectation value of the oscillator coordinate, (q), varies similarly 
to the expectation value of 2{a z ) (the difference between the occupation probabilities of 
the two-level system). In addition, the difference 2(a z ) is stair-case like. This behavior is 
explained by the interplay between the adiabatic region (whenever X(q) » Q) and non- 
adiabatic region (whenever X(q) < Q). The transition from one level to another occurs 
during short periods of time when the system is out of its adiabatic region. Thus, the 
system can be locked on one level if it is prepared to fulfill the adiabatic condition all the 
time. 
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